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ABSTRACT 

We present here the results of the observation of CTB 37A obtained with the X- 
ray Imaging Spectrometer onboard the Suzaku sateUite. The X-ray spectrum of CTB 
37A is well fitted by two components, a single-temperature ionization equilibrium 
component (VMEKAL) with solar abundances, an electron temperature of kT^ ^ 0.6 
keV, absorbing column density of A^h '^ 3 x 10^^ cm~^ and a power-law component 
with photon index of F ^^ 1.6. The X-ray spectrum of CTB 37A is characterized by 
clearly detected K-shell emission lines of Mg, Si, S, and Ar. The plasma with solar 
abundances supports the idea that the X-ray emission originates from the shocked 
interstellar material. The ambient gas density, and age of the remnant are estimated 
to be ~l/~'^/^cm~'^ and ~3x 10^/^' ^ yr, respectively. The center-filling X-ray emission 
surrounded by a shell-like radio structure and other X-ray properties indicate that this 
remnant would be a new member of mixed-morphology supernova remnant class. 

Key v^fords: ISM: supernova remnants-ISM:individual(CTB 37A)-X-rays:ISM 



1 INTRODUCTION 

CTB 37A (also called G348.5-h0.1, RA(2000) = 
17^14'"06°, Dec. (2000) = -38°32') was discovered by 
IClark. Caswell fc Green I (119751 ) in the radio band and has 
a shell-type morphology with an angular size of 15 arcmin. 
From Very Larg e Array observations at wavelen gths of 6, 
20, and 90 cm, iKassim. Baum fc Weilerl (|l99ll ) reported 
that this supernova remnant (SNR) is expanding in an 
inhomogeneous region, and is part of a complex composed 
of three SNRs: CTB 37A, G348.7-h0.3 (also called CTB 
37B), and G348. 5-0.0. Fro m the ASCA Galactic plane 
survey data of G348. 5-1-0.1, lYamauchi et all l|2008l ') found 
that the X-ray spectra of the SNR was heavily absorbed 
by interstellar matter with Ah ~ 2 x 10^^ cm~^ and 
the size o f X-ray em i ssion was comparable to its radio 
structure. iFrail et al.l (|l996l ) detected OH masers with 
velocities at about 20 and 6 km s~^, in the d ir ection 
of CTB 37A at 1720 MHz. iRevnoso fc MangumI (|2000l ) 
surveyed the environment of this remnant with associated 



OH 1720 MHz masers in the CO J=l-0 transition with the 
12 Meter Telescope of the NRAO. They reported that a 
number of m olecular clouds ar e inter acting with the SNR 
shock fronts. lAharonian et al.l ([2003) using Chandra and 
XMM-Newton data showed the presence of thermal X-rays 
from the Northeast part, an extended non-thermal X-ray 
source, CXOU J171419. 8-383023 in the Northwest part 
and a 7-ray source, HESS J1714-385, coincident with the 
remnant. They found a high absorbing column density of 
Ah ~ 3 X lO^'^ cm~'^. They claimed that the observed 
X-ray morphology was a result of interaction with the 
inhomogeneous medium surrounding the remnant and this 
inhomogenei ty was also responsibl e for the break-out radio 
morphology. ICastro fc Stand (|2010|) using observations with 
the Fermi-LAT, have revealed 7-ray emission from this 
SNR; the spectrum of the source coincident with CTB 37A 
was fitted by a power-law (PL) model with an exponential 
cutoff at energy i?cut=4.2 GeV. Considering the lack of 
evidence for contribution by a pulsar and the presence of 
maser emission for the remnant they proposed that 7-rays 
result from SNR-molecular clouds interactions. 
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The distance to the CTB 37A has been estimated 
from 21 cm absorbtion me a surem ent to be in between 6.7- 
13.7 kpc bv ICaswell et al.l (|l975r ). From velocity measure- 
ments of molecular clouds associated with the remnant, 
iRevnoso fc MangumI (|2000l ) adopted a distance of 11.3 kpc. 
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So we will use d=11.3 kpc for our calculations throughout 
this work. 

The location o f CTB 37A containing OH mascr sources 
dFrail et al.1 ll996D . IFGL J1714.5-3830 (Castro fc Sla ne 
I2OIOI ), X-ra y source CXOU J1714 19.8-383023, and HESS 
J1714-385 JAharonian et al.1 l2008t ). as well as two addi- 
tional SNRs (G348. 5-0.0 and G348.7-h0.3) make this SNR 
interesting and important. In addition, its morphology im- 
plies a similarity to the recently proposed group of mixed- 
morphology (MM) SNRs, increasing its importance. High 
quality imaging and spectra obtained from the data pro- 
vided by X-ray observatory Suzaku are used to produce the 
results in this work. 

The organization of this paper is as follows: we describe 
the Suzaku X-ray observation of CTB 37A, including details 
of data reduction in Section 2. The image and spectral anal- 
ysis are given in Sections 3 and 4, respectively. Finally, con- 
sidering its morphology and investigating radial variation of 
the electron temperature, we discuss the implications of MM 
class, the nature of thermal and non-thermal components of 
CTB 37A in Section 5. 



of 0-1.5, 1.5-2.5, 2.5-3.5, 3.5-4.5 arcmin. The extended 
non-thermal X-ray source (CXOU J171419. 8-383023) is ex- 
cluded from spectral analysis as shown in Fig. 1. The 
dashed black circle centered at (RA(2000) = 17*' 14" 19", 
Dec. (2000) = -38°24'36" with radius 1.7 arcmin) repre- 
sents the background region. The lower left corner in the 
FOV that contains calibration source emission is also ex- 
tracted. 

Figure 2 shows the XISO image in 0.3—10 keV energy 
band, whic h is overlaid with the rad io image obtained at 
843 MHz bv lWhiteoak fc GreenI (|l996l ') for comparison. Dark 
crosses indicate the direction of the detected OH (1720 MHz) 
maser emission at velocities ~ —65 km s~^ associated with 
CTB 37A, and the white crosses show the positions of maser 
emission at velocities ~ —22 km s~^ ([Frail et al.ll 19961 ). The 
diamond indicates the position of CXOU J171419.8-383023 
and the circle represen ts the location of HESS J1714-385 
JAharonian et al.ll200^ . 



2 OBSERVATION AND DATA REDUCTION 



Suzaku (JMitsuda et al. ' 2007) is the fifth Japanese X-ray 
astronomy satellite launched on 2005 July 10. Suzaku has 
observed CTB 37A on 2010 February 20. The observation 
ID and exposure time are 504097010 and 5 3.8 ks, respec- 
tively . The X-ray Imaging spectrometer fXIS: lKovama et al.l 
(|2007f )l consists of four sets of X-ray CCD camera system 
(XISO, 1, 2 and 3). XISO, 2, and 3 have front-illuminated (FI) 
sensor and provide coverage over the energy range 0.4—12 
keV, while XISl has a back- illuminated (BI) sensor pro- 
viding greater sensitivity at lower energies (0.2—12 keV). 
The XIS has a field of view (FOV) of 17.8x17.8 arcmin^ 
(1024 X 1024 pixels). Each XIS CCD has an ^'^Fe calibration 
source, which can be used to calibrate the gain and test the 
spectral resolution of data taken using this instrument. The 
XIS2 sensor was available only until 2006, therefore, we use 
dataofXISO, XISl, XIS3. 

Reduction and analysis of the data were performed fol- 
lowing the standard procedure using the headas software 
package of version 6.5 and spectral fitting was performed 
with XSPEC version 11.3.2 (Arnaud 1996 ). The XIS was op- 
erated in the normal full-frame clocking mode, with the stan- 
dard 3x3 and 5x5 editing mode. We generated XIS re- 
sponse matrices using the xisrmfgen software, which takes 
into account the time-variation of the energy response. As 
for generati ng ancillary response files (ARFs), we used xis- 
SIMARFGEN (|lshisaki et al.l 120071 1. The latest version of the 
relevant Suzaku CALDB files were also used. 



3 IMAGE ANALYSIS 

Figure 1 shows XISl image in 0.3—10 keV energy band. 
We extracted the spectrum from the brightest region 
represented with the outermost solid circle centered at 
RA(2000) = 17''14"'30', Dec. (2000) = -38°32'07" with 
a radius of 5.5 arcmin. To derive the radial variation of the 
electron temperature kTc, we take four apertures with sizes 



4 SPECTRAL ANALYSIS 

XIS spectra were extracted using xselect version 2.4a 
from all the XISs with a circular extraction region of radius 
5.5 arcmin and are grouped with a minimum of 50 counts 
bin~^. We fit the spectra with a coUisional ionization 
equilibrium (CIE) mo del with variable abundances (xspec 
model 
(1981); 
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JLiedahl. Osterheld fc GoldsteinI ( 19951)) modified by inter- 
stellar absorbtion (wabs in xspec, iMorrison fc McCammonI 
(|l983r )). The parameters of the absorbing column density 
(A'h) and electron temperature (fcTe) ar e set free while all el- 
emen ts were fixed at solar abundances IjAnders fc Grevessd 
[l989i). The best-fitting reduced x^/d.o.f. for this model is 
2222.5/832= 2.67. To find out if there was any contribution 
from non-thermal emission we added a PL component 
(VMEKAL-hPL) yielding a better reduced x^ of 935.5/830 
= 1.13. Then, Mg, Si, S and Ar lines in the spectrum were 
set free, while the rest were fixed at their solar values, we 
found insignificant improvement in y^ value (891.6/826 — 
1.08). Therefore, we decided to fix the abundances of Mg, 
Si, S, and Ar at their solar values. In Table 1, we present 
the best-fitting parameters and the statistics obtained with 
an absorbed VMEKAL-I-PL with corresponding errors at 90 
per cent confidence level (2.7 o). Figure 3 shows the spectra 
of the XISO, XISl, and XIS3 simultaneously, in the energy 
range of 0.3—10 keV that is taken from the region shown 
by the solid dark circle (the outermost) presented in Fig. 1. 
On the other hand, we performed annular spectral analysis 
for four regions shown by the circles in Fig. 1 to be able to 
derive the radial temperature variations of CTB 37A. The 
annular regions are spaced by 1 arcmin from the innermost 
circle with radius r=1.5 arcmin. VMEKAL-I-PL spectral 
model were also fitted to each annulus, while fitting them 
we kept the absorbing column density A^h at its best-fitting 
value for the entire remnant. Figure 4 shows the electron 
temperature variations with respect to the radius. 
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Table 1. Best-fitting parameters and x^ values of the spectral fitting in the full energy band (0.3—10 keV) 
with an absorbed VMEKAL+PL model with corresponding errors at 90 per cent confidence level (2.7 a). 



Component 



Parameters 



Wabs 

VMEKAL 

Abundance" 



NuixW^ 



') 



PL 



fcTc(keV) 

Mg 

Si 

S 

Ar 

VEM^ 

Flux£" 

Photon Index 

norm(xlO~^photons cm~^s 

Flux_£ 



VMEKAL+PL 


2.9 ±0.1 


0.63 ±0.02 


(1) 


(1) 


(1) 


(1) 


72.3 ±2.1 


1.36 ±0.03 


1.6 ±0.1 


1.9 ±0.2 


1.5 ±0.1 



935.6/830=1.13 



GrevessdlToSgh . 



"^ (1) indicates that the elemental abundance is fixed at solar 11 Anders &: 

Volume emission measure VEM= / nenndV in the unit of 10®* cm~^, where He 
densities of electrons and protons, respectively, and V is the X-ray emitting volume. 
"^ Unabsorbed flux in the 0.3 — 10 keV energy band in the unit of 10~^ erg s~^cm~^, 
"^ Total unabsorbed flux of the sum of the VMEKAL and PL components in the 0.3 - 
in the unit of 10~^ erg s~^cm~-^. 



and rtH are number 



10 keV energy band 



5 DISCUSSION AND CONCLUSIONS 

In this work, we provide a description of the X-ray emission 
of CTB 37A based on Suzaku archival data. We obtained a 
clear image and high quality spectra of diffuse X-ray emis- 
sion. We have examined the thermal and non-thermal emis- 
sions coming from the remnant. The X-ray spectrum of CTB 
37A is characterized by thermal emission dominated by K- 
shell emission lines of Mg, Si, S, and Ar which are clearly 
detected. 

As seen from Fig. 2, the radio emission of CTB 37A 
comprises a partial shell towards the north and east and 
an extended outbreak to the south, while the X-ray emis- 
sion has a deformation along the southwest limb, where the 
morphology appears indented. The reason for such a defor- 
mation may be due to the inhomogeneous medium along 
this specific region and this may well be supported by the 
fact that remnant is close to the Galactic plane and several 
OH masers at 1 720 MHz are detected towards CTB 37A 
l|Frail et al.lll996l ) (shown by crosses). It may also be relevant 
that 7-ray emission thought to be associated with the inter- 
action between CTB 37A and dense surro unding material 
that has been detected with the Fermi LAT (jCastro fc Sland 
I2OIOI ). 

5.1 Implications for mixed-morphology 

SNRs were originally divided into shell-like. Crab-like (ple- 
rionic) , and compos ite (shell- like containing plerions) rem- 
nants (|Sewardlll985t ) according to their X-ray morphology. 
Recently, an additional MM class (also called thermal com- 
posite) appeared, which are center-fill ed in X-rays and shell- 
like at radio wavel engths |Seward 19901 : I Jones et al.l 119981 : 
iRhofc Petr'3ll998l ). As seen from Fig. 2, CTB 37A has 
a shell-like morphology in the radio band while centrally- 
filled in X-ray band, in this regard, our first impression is 
CTB 37A seems to be a MM SNR. Examples of well-known 



MM SNRs include W 28 (iRho fc Borkowskill2002l). G290.1 - 
0.8 jS lanc ct al."2002'), arid IC 443 (JKawasaki etal.1 120021 '). 
[Lazendic & Slanc (2006) reported important results about 
this class by compiling a list of 26 MM SNRs. 

The X -ray char acteristics of MM remnants have been 
defined by iRho fc P ctrc (1998) as foUows: (1) The radial 
temperature distribution is relatively flat; (2) The X-ray 
emission arises primarily from shocked interstellar mate- 
rial, and not from ejecta; (3) The remnants are typically 
located close to molecular clouds or very dense regions; 
and (4) The dominant X-ray emission is thermal in na- 
ture. Subsequent studies indicated that MM SNRs had a 
complex plasm a structure with multiple components (e.g. 
Rho fc Borkows ki (200 2.)) and enhanced abundances (e.g. 



Yamauchi et al.l (|1999l 'l: ISlane et al.1 (|2002l ')'). and they have 
evolved over ^10* yr, which means that the plasma is i n 
CIE or an overionization condition (jKawasaki et al.ll2005l ). 
As can be seen from Fig. 4, the annular analysis of CTB 
37A may well be indicating a very small scale radial vari- 
ation in its temperature between the selected regions. The 
best-fitting metal abundances are found to be solar in gen- 
eral, confirming the absence of ejecta contamination in se- 
lected regions. It may support the idea that the X-ray emis- 
sion originates from the shocked interstellar material. The 
plasma of CTB 37A is in a coUisional ionization equilib- 
rium condition and is located in a region with density varia- 
tion, possibly associated with molecular clouds. The plasma 
has thermal and non-thermal emission, but the emission is 
dominated by thermal component (~ 90 per cent of the 
total X-ray flux). These X-ray properties of CTB 37A ex- 
emplify the typical characteristic of MM SNRs as defined by 
iRhofcPetrdl 19981 1. 

There are a few models that can produce centrally en- 
hanced X-ray emission such as evaporation of clouds left 
relat ively unspoiled a fter t he passage of the SNR blast wave 
(e.g. IWhite fc Longj (|l99ll )) and "fossil" thermal radiation 
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that is detectable as thermal X-rays from the hotter interior 
as the shell of an expanding SNR cools below ^10 ^ K and 
becomes invisible due to interstellar absorbtion (e.g. ISeward' 
l|l985ri l. as the SNR evolves, the temperature and density of 
the hot interior plasma gradually b e come uniform through 
thermal conduction (e.g. ICox et aU ([1993)) or evolution in 
a med ium with a d ensity gradient viewed along the line of 
sight (|Petrukll200lh . The evaporation model requires dense 
clouds, the thermal conduction model requires a relatively 
high density ambient medium. CTB 37A is located in a re- 
gion of greatly varying density, with OH maser sources in- 
dicating interaction with molecular clouds. In this regard, 
the center-filled X-ray morphology of CTB 37A is consis- 
tent with evaporating clouds model. The radial temperature 
variation in the plasma of CTB 37A (fcTc ~ 0.6- 0.8 keV) is 
consistent with o ther MM SN Rs such as W44 fRho et al.l 
1994; Shclton. K untz fc Petrel [2004). 3C391 (Rho & Pctd 
19961 : IChen et all 120041 ') and HB21 (|Pannuti et al.l 1201(1 )7 



Very small temperature variation in the plasma of CTB 37A 
as shown Fig. 4 can be explained by both evaporation and 
thermal conduction models. Future deep X-ray observations 
and detailed spectral analysis of this remnant would give 
more detailed information to compare with theoretical mod- 
els that produce MM SNRs. 

5.2 Thermal component 

The X-ray emission of CTB 37A is dominated by thermal 
emission that can be best described by an absorbed CIE 
plasma model (VMEKAL) with an absorbing column den- 
sity of A^H ~ 3 X 10^^ cm~^, an electron temperature of 
kTc ~ 0.6 keV, and solar abundances of Mg, Si, S, and Ar, 
which indicate a shocked interstellar/circumstellar material 
origin. 

For full ionization equilibrium, the ionization timescale, 
r = rici, is required to be > 10^^ cm~^s, where t is the 

J ilasma age or the time since the gas was shock-heated 
Mas ail [l98J). To determine the age of the remnant, ric 
should be estimated from the emission measure, ncnuV , 
which is related to the normalization of the VMEKAL model 
according to the equation, norm=nonHV/(47rd'^10^^), where 
V is the X-ray emitting volume, tih is the volume density of 
hydrogen and d is the distance. For simplicity, we assumed 
the emitting region to be a sphere of radius 5.5 arcmin. 
Considering the possibility that less than the entire volume 
is filled, we write the volume V=Ki./, where Vs is the full 
spherical volume, / is the filling factor. We then carry the / 
factor through our calculations to show the explicit depen- 
dence of each derived quantity on this factor. Knowing that 
the SNR is at a distance of 11.3 kpc and Ue = 1.2nH, we 
estimated the emission volume to be F ~ 6.7 x 10^^/ cm''. 
Consequently, we find an ambient gas density of ~l./~^ 
cm~^ and age of ~3 x 10*/^'^^ yr (assuming Uct ~ 1 x 10^^ 
cm"^s), implying that CTB 37A is a middle-aged SNR. Fi- 
nally, we calculated total mass of the X-ray emitting plasma, 
Mx, by Mx=7nHWcl^~ 53O/^"M0, where jtih is mass of a 
hydrogen atom, /i=0.604 is the mean atomic weight. 



5.3 Power-law component 

The Suzaku X-ray spectral data of CTB 37A is well fitted 
with a thermal component and an additional hard compo- 



nent. There could be a few reasons for the hard X-ray emi- 
sion: (i) an association with a classical young pulsar, (ii) 
a contribution from an extended non-thermal X-ray source 
(CXOU J171419.8-383023), (iii) overionization of the plasma 
which pro duces excess hard emiss ion as has been the case 
for IC443 (|Yamaguchi et al.l 120091 ). The hard component is 
well fitted by a PL model with a photon index value of '^ 1.6. 
This value is consistent with that of cla ssical young pulsar 
value ranging in between 1.1 and 1.7 (|Chakrabartv et al.l 
l2001f ). However, there is no pulsar reported that is asso- 
ciated with this remnant. In the Northwest region of the 
CTB 37A, an extended non-thermal X-ray source (CXOU 
J171419. 8-383023) is reported (RA(2000) = 17''14'" 20°, 
Dec. (2000) = -38°30'20") bv lAharonian et al.l l|2008l ). In 
their work, a non-thermal emission from the source with a 
spectral index of ~1.32 is found which is lower (harder) than 
our best-fitting value of ~1.6. Although we excluded CXOU 
J171419. 8-383023 (with radius 2.1 arcmin) from our spectra 
during our spectral analysis, our fits required a non-thermal 
component. To investigate this, we performed spectral anal- 
ysis also for individual regions by selecting small rectangu- 
lar regions that are being further away from the known ex- 
tended non-thermal source. The non-thermal flux is found to 
be stronger for the selected small regions nearby the source 
compared to the ones further away. We have obtained an 
unabsorbed flux value of F^ ~ 1.4 x 10~^ erg s~^cm~^ for 
the non-thermal extended source in the 0.3—10 keV energy 
range. When we compare it with the unabsorbed flux value 
(Fx ~ 0.14 X 10~® erg s~^cm~^) of PL component of best- 
fitting, we find a factor of ten difference between them. This 
difference indicates that the extended source is most likely 
the origin of the PL component and emission scattered from 
the source into the field of the rest of the remnant by a broad 
point spread function of the Suzaku mirrors. 

The spectral studies of CTB 37A indicate that the 
plasma is best described by a thermal component in CIE 
condition with solar elemental abundances and a non- 
thermal component with a photon index of ~1.6. Thermal 
emission possibly originates from the shocked interstellar 
material with ambient gas density of ~l/~'^"cm~''. The 
best spectral fits require an ionization timescale of r ^ 10^^ 
cm~''s, implying an age of ~3xl0*/^'^ yr. The origin of the 
power-law component is more likely the effect of the contri- 
bution from the extended source (CXOU J171419.8-383023) 
located in the Northwest part of the remnant. CTB 37A is 
most likely a new member of mixed-morphology SNR. 
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R.A. 

Figure 1. Suzaku XISl image of CTB 37A in the 0.3—10 keV energy band. The spectral integration regions for the source (annular 
apertures) and the background are indicated with the solid dark and dashed black circles, respectively. The X-ray source (CXOU 
J171419. 8-383023) is extracted from the spectral analysis as shown in the figure. The corner of the CCD chip illuminated by ^^Fe 
calibration source is excluded from the image. The coordinates (RA and Dec) are referred to epoch J2000. 
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Figure 2. CTB 37A is displayed in 0.3—10 keV energy band X-ray data obtained from Suzaku XISO detector. The image is overlaid 
by 843 MHz radio isocontours logarithmically scaled (0.23, 0.36, 0.57, 0.90, 1.42, 2.23 and 3.51 mJy/Beam). Radio emission clearly 
shows 2 radio SNRs, CTB 37A and G348. 5-0.0. X-ray data is smoothed with 3 pixel Gaussian in order to highlight the structure. The 
colour-coding shows the brightness levels from 1.3 counts/pix to 13.36 counts/pix in logarithmic scale from bright to dark. The masers 
from the field are indicated with dark and white crosses at velocities ~ —65 km s~^ and ~ —22 km s~^, respectively. The diamond 
shows the extended non-thermal X-ray source CXOUJ 171419.8-383023. The location of HESS J1714-385 is represented with the circle, 
in which the positional error is specified with the radius of 1.2 arcmin. The coordinates (RA and Dec) are referred to epoch J2000. 
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Figure 3. Background-subtracted XIS (XISl, XISO and XIS3) spectra of CTB 37A in 0.3-10 keV energy band fitted witii an absorbed 
VMEKAL and PL model. Tlie bottom panel is the residuals from tlie best-fitting model. 
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Figure 4. The radial variation of the electron temperature of CTB 37A. 



